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Abstract 

The possibility is discussed of using straight-line paths of integration in computing the integral 
representation of the three-body Coulomb Green's function. In our numerical examples two dif- 
ferent integration contours are considered. It is demonstrated that only one of these straight-line 
paths provides that the integral representation is valid. 

PACS numbers: 03.65.Nk 



*E-mail: 



zaytsev@fizika.khstu.ru 



I. INTRODUCTION 

In two previous papers the method was introduced as a new approach for solution 
of the three-body continuum problem using infinite set of L 2 parabolic Sturmian basis func- 
tions for the wave function of the system. The goal of these papers has been the construction 
of exact analytic matrix elements of the three-body Coulomb Green's function. The cor- 
responding six-dimensional resolvent operator has been expressed as a convolution integral 
of three two-dimensional Green's function. In this paper we wish to learn how to choose 
appropriate straight-line paths of integration which provide that the integral representation 
is valid. 

Below we outline how the Schrodinger equation for a three-body Coulomb system is trans- 
formed into a Lippmann-Schwinger equation in terms of generalized parabolic coordinates. 

The Schrodinger equation for three particles with masses m 1; m 1; m 3 and charges Z 1; Z 2 , 
Z 3 is 
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where R and r are the Jacobi vectors 
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removes the eigenenergy E = ^— K 2 + k 2 giving the equation for \I> 
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Then, the operator in the square braces is expressed in terms of the generalized parabolic 
coordinates 3] 
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where k; s = k ' ms , kam ' is the relative momentum, k; s = k ^-, h H = Ik; J. In the resulting 

'* mi+m s ' 18 ki s ' I I ° 

equation 
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the first operator is given by 
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D^Yl 77" V h^ + h Vj +2k ls ti s , (8) 

for j ^ I, s and I < s. Here t/ s = Zl ^ sftu ; = ; the one-dimensional operators ft,^ 

and /i^ are 

d ^ d d \ t „ f d d ., 9 



Do is the leading term which provides a three-body continuum wave function that satisfies 
exact asymptotic boundary conditions for Coulomb systems, when the three particles are far 
away from each other j^]. In turn, the operator D\ (which contains the non-orthogonal part 
of the kinetic energy operator) is regarded as a small perturbation which does not violate 
the boundary conditions. 

The best known approximate solution to the equation (HQ), the so-called C3 model 13, 

nnn * 

m [5|, |6j, is obtained by neglecting of D\. Many improvements to the C3 model have been 
developed by considering in some approximate way of the neglected terms of the kinetic 
energy (see, e. g., [7j and references therein). In our approach the wave function \I/ is 
obtained by solving an equivalent Lippman-Schwinger integral equation. We multiply (JTj) 

3 

by EI Ms (6j + Vj) from the left before the transformation of @ into an integral equation. 

3=1 

Thus, in the resulting equation 

iF = !7 {0) - fevy (10) 

& plays the role of Green's function operator which is formally inverse to the six- dimensional 
operator \) given by 

3 

h = IT P-U (& + Vj) D = A*13 (6 + Tfy) fit2 (6 + Vs) f)l 

(11) 

+P23 (fl + Vl) A*12 (& + ffe) $2 + ^23 (fl + Vl) A*13 (6 + Tfe) ^3, 

= ^. + h v . + 2k ls t ls . (12) 
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The "potential" V is denned as 



(13) 



The inhomogeneous term W ' of Eq. (JTDD can be taken as the wave function of the C3 
model, i. e. expressed in terms of a product of three Coulomb waves. 

It has been suggested in to treat the equation within the context of L 2 parabolic 

Sturmian basis set 



i=i 

<f>n jmj (&, vj) = fa, (s,) '■;„: (///) 



(14) 
(15) 

%l) n (x) = V2be' bx L n (2bx), (16) 
where b is the scaling parameter. A solution \I> of the Lippman-Schwinger equation (|10p is 



expanded in basis ( 1141) as 



(17) 
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The discrete analog of the Lippman-Schwinger equation is obtained by putting (fTOj) in the 
basis set f fT4|) . This yields 

a = a (o) -0_Ka, (18) 

where and V are the operators G5 and V matrix representations in basis (114j) . a and g^ -* 
are the coefficient vectors of \I/ and respectively. 

It has been shown in our previous paper ^ that the matrix (5 can be represented in the 
form of a convolution integral 



<* (±) = v&?J I fig G (±) (** ® GW (* 13 ; £ 2 ) 

C(!)C(2) 
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Here (ij s , £j) is the matrix which is inverse of the two-dimensional operator 



(20) 



fy? + ( — ) (0 + ma tri x representation in the basis (fT5|) . i. e. 
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& Q, 



G (±) (t ls , £j)=lj. 
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Here the matrix hj of the operator fyj (fl2|) is expressed in terms of the one-dimensional 
operators Qj matrices: 



hj = h ?j <g> I Vj + <g> h % + 2/^^-. (21) 

In Eqs. f )20l) and fl2Tl) I*. , 1^. and I,- = I £ . are the unit matrices. Qj = ®I Vj +I& ®Qjj, , 
where Q^. and Q^. are the matrices of ^ and r]j in basis (TToT) . respectively. 

In this paper we make use of matrices (£ £) flA"20l) which are more symmetric (in £ 
and 77) than that obtained in our previous work [2j. For simplicity of the notation, we omit 
indices for a while. The new matrix, e. g., G^ (t, 8) also obeys the completeness relation 

/ d£G (+) (t; £) = [Q £ ®I^ + ^® Q v r l , (22) 

Z7TZ J 
C 

established in 2|. Here C is a contour originating at S = 00, below the positive real 
axis rounding the lowest bound state Si = for t < (or the origin for t > 0), 

and then heading back to £ = 00 — this time staying above the cut (see Fig. 1). The 
integration contours of the convolution integral (j!9p C^ 1 ' 2 ^ are similar to the contour C 



(see e. g. {9]). However, despite the known paths of integration the representation (fl9|) 
poses several difficulties in practical applications, the most serious of which is that one 
cannot trace crossing the cut along the positive real axis in ^3-plane during integration 
over £ 12 . The (numerical) evaluation of the integral ffl9l) can be simplified considerably by 
using straight lines as paths of integration. In this paper we wish to learn how to choose 
appropriate straight-line paths and for which the integral representation (TTTJ]) is 
valid. Unfortunately, the contour integrals of interest cannot be treated analytically, so we 
must resort to numerical experiments. 

The numerical examples presented in Section III show that the contour C in (122]) could be 
deformed so that it becomes the disconnected pair of straight lines. The value of the integral 
over each of these straight-line paths is half the value of the contour integral (l2"2l) . In this 
section we consider two kinds of straight-line paths. In Section IV based upon the numerical 
results obtained for double integrals which arise from the matrix product we find 

the straight-line contours C^ 1,2 ) providing a non-zero integral representation ffl9l) . Section 
V contains a brief discussion of the overall results. For completeness we review briefly the 
results of our previous works [l, ^ in the Appendix. 
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II. PRELIMINARIES 



We assume that the relationships between contour integrals obtained in this work re- 
sult from the Green's functions properties and are independent of the base function (|14p 
numbers (that specify indices of the matrix elements of the operators). Thus in all our 
numerical examples (except for the case where the inverse relationship between f) and (3 
is demonstrated) we only use the element Gq o- 0,0 falsi £j) = @o falsi £3) °f the matrices 
G( + ) (ti s ; Sj). Notice that the completeness relation (1221) for Gq (t; S) takes the form 

^- J d£G { + \t; 5) = 26. (23) 



c 



We consider two electrons in the Coulomb field of He ++ , i. e. Z\ = Z 2 = — 1 and Z 3 = 2 
(we use atomic units hereafter). Infinite mass for the nucleus is adopted. Let the two 
electrons move in opposite directions with equal energies E\ = E2, i. e. ki3 = k and 
k 23 = — k. Note that in this case k 12 = | (ki3 — k 2 3) = k (so that matrix elements of the 
"potential" operator V ([TBI could be evaluated analytically). Further we set the electron 
energies E\ = E2 = 25/Ry and the scaling parameter 6=1. 

III. STRAIGHT-LINE PATHS 

Clearly, the contour C can be deformed by moving its left edge toward minus infinity so 
that it becomes the disconnected pair of straight lines (Ci and C 2 in Fig. 1) parallel to the 
real axis. This deformation does not alter the value of the integral 



m: m! m! 



J n m; n' nv 

27TI 

C 



2^ / d ZG { :l n , m ,{t- £). (24) 



Hence v nm;n ' m ' is divided into two contributions: 

Vnm;n'm' = v n m; ri m' V nm; ri m' i (^5) 

where 



nm;n>m> 2 m 



l - j d£G { n + i, n i m >(t; £), (26) 

Ci 

C 2 



Consider the integral (I26I) . Notice that the energy 8 is parametrized by 8 = iyo + x on 
the contour C\ . Obviously, the value of v^' m . n , m , is independent of the parameter y which 
specifies the position of C\ with respect to the real axis (see Fig. 1). The numerical results 
for the integrals 



oo 



V ° ] = J dX ^n Gi ° +) \P W ° + X ) (28) 

— oo 

with different yo are presented in Table 1. Note that the integrand in (|28|) involves hypergeo- 
metric functions 2F\(a, b; a + b; z). We found that the Gauss continued fraction 10J (rather 
than the infinite sum (15.3.10) in |ll|) provides an efficient approximation for the hyper- 
geometric functions. The integrals are computed using IMSL FORTRAN Library routines. 
The real and imaginary parts of the integrand identified as A[x) and B(x) are displayed 
in Fig. 2 for different y^. It is seen in Fig. 2 that B(x) tends to zero as the parameter y 
increases uniformly with respect to x. This property of B(x) is consistent with the negligible 
imaginary parts of the integrals in Table 1. Using this observation, and (1231) . we conclude 
that 

=vf ) = l = \v Q . (29) 
Moreover, our numerical computations show that 

V n m; ri m' = V n m; n' m' = 7} Vn m '< n> m ' 

holds for arbitrary n, m, n', m'. 

Another straight-line path is obtained by rotating the contour C\ about some point Xq on 



the positive real axis through an angle (p in the range — 7r < (p < [2|, [l2| . For definiteness, 



we choose xq = | and ip = — ~, i. e. 8 = | + iy on the resulting contour C3, shown in Fig. 3. 
The path C 3 crosses the cut so that its lower part (depicted in Fig. 3 by the dashed line) 
descends into the "unphysical" sheet (— 2tt < arg(£) < 0). To analytically continue matrix 
elements of g^- 1 and g^+J in ( 1A20I) onto the unphysical sheet we use the formulae ( 1A16I) 
and flA17l) . Note that the numerical result for the integral 



-00 



4 i] = l [dyG^(-l- ^+iy) (31) 
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presented in Table 1 also satisfies 



i 3) = lv . (32) 



Thus, we find another straight-line path C3 for which 



,(3) 



n m; n m' 



(33) 
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IV. DOUBLE INTEGRALS 

Now that we have the relationships (I3"U1) and fl3"31 between the integrals along the contour 
C and the integrals over the straight-line paths, we can use C\ and C3 in the integral repre- 
sentation (|T9|) of the Green's function operator. Before proceeding, we consider the matrix 
product f)0 ( - + ' 1 to determine the normalizing factors K corresponding to the paths C\ and 
C3. Using the six-dimensional operator (3 (fill) matrix representation 



f) = /Ul3AH2hi ® Q 2 ® Q 3 + /U23AH2Q1 ® h 2 <g> Q 3 + A*23^13Ql ® Q2 ® h. 



3) 



(34) 



(HHD and PJ, we have 



where 



fj = K {Wi + W 2 + W 3 } 



Wi 



1 Ml 2 

(27Tl) 2 /^23 



/ / dS l dE 2 Ii ® [Q 2 G« (t 13 ; £ a )] 

C« C( 2 ) 

b 2 . / b 2 \ . / b 2 



Q 3 G(+) (t 12 ; ^ + g I** - £ 1 ) + -E, 



M13 V 2 



Wo 



1 M12 
(27ri) 2 /ii3 



/ / d£ lC f£ 2 [QiG«(t 23 ; £1)] ®I 2 



CC 1 ) c« 

Q:sG , + ) (ti2 . M, + g _ £l ) + Mia (Ma _ s 2 



M13 V 2 



W 3 = — K- J J d£ x d£ 2 [QiG {+) (t 23 ; S x )\ ® [Q 2 G« (fo; £2)] ® I 3 



(35) 



(36) 



(37) 



(38) 



The matrices f) and (S^-* must be inverses of each other. Therefore, the normalizing factor 



K and the matrices W,-, j — 1, 3 satisfy the condition 



N {Wi + W 2 + W 3 } = I. 



(39) 



If we choose C^> = C x and = C x (or = = C 3 ), then it follows from fl22} and fl3U]) 
(or d33])) that 

W 3 = Jli®I 2 ®l3= (40) 



Further, we assume that 



w 1 = -i, W 2 = |I, 



(41) 



so that the sum {Wi + W2 + W3} is proportional to the unit matrix I. In turn, the 
constants a and (3 can be determined by, e. g., the ratios 



where 



a = — , P = —, 42 
w 3 w 3 



W ^ = WS I I d£ 1 d£ 2 Gi +) (t 13 ; £ 2 ) 
xG« (t 12 ; f + g (f - fc) + ^ (f - £ 2 ) 

= (WE / / df^GW (t 23 ; £1) (43) 

y^W f+ ■ fc ?2 I M12 ^23 _ C \ , M12 A13 _ C 

^3 = 7^ / / ^1^2 G { +) (t 23 ; f 1) 4 +) (tia; £2) 



(2tu) 2 

C(i) C( 2 ) 

Thereafter, the normalizing factor K is expressed as 



N = -. (44) 



First we consider the path C\. 



a) C« = C( 2 ) = d 

In this case the energies £\ and £ 2 are parametrized by E\ = iyo + x\ and £ 2 = iyo + x 2 
with y = 100. Assuming that the energy £ 3 = k 2 — zy — | (#1 + #2) lies 011 the "physical" 
sheet i. e. < arg(£ 3 ) < 2tt, in view of (1251) - ([291 . we obtain that 



00 , . 00 



h I dxG?> (i; k 2 -^y Q -\x)=2±- ^ J dxG +> (£; tf-iifo + s) 

(45) 

c 2 

Hence, one might expect that the constants a and /3 ( 1421) are negative. Note that in our 
case the integrals Wj ( T43j) take the forms 

Wl = ~J^f\ 1 1 dx idx 2 G +) (-§ ; %o + s 2 ) Go +) (i; /c 2 -%yo~\ {x x + x 2 )) , 

OO 

w 2 = -j^l J J dxidx 2 G { +) (-§; ij/ + x x ) G { +) (i; A; 2 -iy -\ (x 1 + x 2 )) , (46) 



(2^)2 



w 3 = -7^)2 / / cfeidx 2 (-|; iy + xi) Gq +) (-|; iy + x 2 ) . 



-00 
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From the numerical results for the double integrals (1461) presented in Table 1, it follows that 
wi = W2 = — \ and u> 3 = 1, i. e. a = (3 = — |. Thus, we have a + (3 + 1 = 0, and 
so the equation ( 1441) is meaningless. This outcome is consistent with the numerical result 



J 0,0 



obtained for the integral in the expression ( fl9i) for the diagonal matrix element of 
corresponding to the basis function |0) = \rij = rrij = 0, j = 1, 3) ( 1T4"1) : 

oo 

J o = ?2^F / / ^1^2 (-f ; iy + x x ) G^ (-§; iy + ^a) 

— ' (47) 



x 



G o (A; fc2 - *yo - 1 (a?i + ^2)) 



which is presented in Table 1. Therefore, we conclude that the contour C\ cannot be used 
in the integral representation (fl9l) . 
Now consider the the contour C 3 . 



b) C« = c( 2 ) = C 3 

The energies £1 and £2 are given by E\ = -y + zyi and £2 = % + iy% on the contour C 3 . 
In turn, the energy £ 3 is parametrized as £3 = -y — f (2/1 + 2/2)- I n this case the contour 
integrals (1431) are transformed into the double integrals 

w 1 = ^ 7 fjJdy 1 dy2Gi +) (-£; £ + zy 2 ) G +) (i; £ - | ( yi + y 2 )) , 

+00 

» 2 = ^f//*AGi +) (-!i T + »fi) ^o +) T - I (fi + W)) . ( 48 ) 

+00 V / \ / 

= j^i J J d yi dy 2 (-f; f + G« (-1 ; f + iy 2 ) . 

+00 

Notice that from (|3T|) it follows that 

—00 

+00 

Therefore, in contrast to the previous case, W\ (w^) would be expected to have the same 
sign as W3. From the result for the numerical evaluations of the integrals (flBl) . presented in 
Table 1., it follows that W\ = w 2 = § and w 3 = 1. Hence a = (3 = | and K = 1. 

To verify that the straight-line path C 3 does provide the desired result, we evaluate 



1 

2^ 



dyGi +) 



1 

2k' 



+ iy 



2v, 



(3) 



(49) 
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numerically the matrix elements 



J 9t, 9T 

xG {+) , , 

«2 "l2i % Win 



(50) 



2. fc^ 
'fc> 2 



and calculate the matrix product t)<3^ of finite size. Notice that the matrix {) is 
"tridiagonal" , i. e. for each pair of indices yrij, n^j and {nij, m'jj, j = 1, 3, the elements 



vanish unless n,- 



< 1 and |m.j 



< 1. Therefore, the minimal rank M m i T . 



of the matrices f) and 0., with which the relation f)(3 = I could be verified, is given by 
Nmin = 2 6 = 64. Actually, to test this equality, we must use all the basis functions |9t) (IT4"|) 
with each of the rij and taking the value one or zero. Our prime interest here is with 



the values of the first row elements 



(+) 



J 0,01 



in the matrix The numerical result 



for the first diagonal element 



(+) 



J 0,0 



, presented in Table 1, corresponds to the inverse 



relationship between () and (5. In turn, the remaining (zero) elements of the first row are 
found to be of the order of 10~ 7 . 



V. CONCLUSION 



In this paper we focus attention on the three-body Coulomb Green's function opera- 
tor representation. The development of our method is based primarily upon the fact that 
for large particle separation the Schrodinger equation is separable in terms of generalized 
parabolic coordinates. Thus, the corresponding six- dimensional resolvent operator can be 
expressed as a convolution of three two-dimensional Greens function. This representation 
includes integration along contours, which encircle the spectra of two-dimensional wave op- 
erators. Unfortunately, these (double) contour integrals are very inconvenient for numerical 
computation. Clearly, it is preferable to employ straight-lane paths of integration. In this 
paper we demonstrate numerically, with two simple examples, that use of an appropriate 
straight-line path of integration provides a non-zero integral representation of the Green's 
function operator. 
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Appendix A: MATRIX REPRESENTATIONS OF ONE- AND TWO- 
DIMENSIONAL OPERATORS 



1. One-dimensional operators 

The matrix representation of the operator f)£ ([9]) in the basis {^(Ol^o (©' which 
is orthonormal with respect £ G [0, oo), is tridiagonal with nonzero elements 

hi in = b + ik + 2bn, hi in _ 1 = (b-ik)n, hi >n+1 = (b + ik)(n + 1). (Al) 

In addition, the symmetric matrix of the operator £ in the basis (fT6l) is also tridiagonal: 



Qn, n' 



Hence, the one-dimensional operator 



_n_ l _ i 

2b' '* ± > 



2n+l 
2b ' 

n+l „l 



n ' — n, 



(A2) 



^, n' = n + l. 



f)^ + 2kt + /xC£ also has the tridiagonal matrix rep- 
resentation [h^ + 2ktl^ + /iCQg] (Ig is the unit matrix) in the basis set (1161) . Fortunately, 
the inverse of the matrix [hg + 2ktl^ + /iCQ^] can be obtained analytically. The elements 
of the resulting matrix g^ ± - ) are expressed in terms of well-known special functions 



7 S(+) 

'711,712 
7«H 







-n 


Qnx-712 














- n 










(ri2 



where n< is the lesser of n\ and n2, and n> the greater of the two. Here 





2b + i(-f-k) 2b-i(-f + k) 

A — — - — , C 



2b-i(-y-k)' 



2b + 1(7 + fc) 



A 

0' 



A; / i 



(A3) 
(A4) 

(A5) 
(A6) 
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and 



k 2 i 2 

^C = --£, £ = -L, (A7) 

(—l) n T in + - — ir) ( 1 1 \ 

IhAj: c, ) - y —l f7lj— j 2Fl ( v ~ n ' 2 + ' T; ~ n + 2 + C J (A8) 



ni 

are the polynomials of degree n in r which are orthogonal 



C m V c 



with respect to the weight function, 



drp(r; ()p n (r; Op m (r; () = 5 nm , (A9) 



^ c) = i r (5 + !T ) r G- !T ) ( - c),T+i (A10) 

(±) 



In (lAlOj) it is considered that | arg(— £)| < n. The function <]4 is given by 



qtHr; C) = (- 



I)" ,' A 2*1 -±»r,n+l; n + - ± *r; C Tl . (All) 



Similarly, the nonzero elements of the tridiagonal matrix of the operator f} v ([9]) in the 
basis {ip n (v)}n=o CED are 

hl n = b-ik + 2bn, h\ n _ x = (b + ik)n, hl n+1 = {b - ik){n + 1). (A12) 

Obviously, the nonzero elements of the matrix of the operator 77 in the basis (FTO]) 
are defined by (IA2I) . Thus, the elements of the matrix g 7? ^ ± - ) , which is inverse of 
the matrix [h^ + 2kt\ v + fiCQrj] (I v is the unit matrix) of the one-dimensional operator 



f) v + 2kt + nCr] 



are expressed as 



gtl(T\ 7) = ^ f^J ^-J^to 0?£V; 0, (M3) 

^7i(r; 7) = ^ ^7-Pn<(r; C) C >+1 C)- (A14) 

Note that (qk ~M are defined in the upper (lower) half of the complex 7 plane where 
ICI > 1 (|C| < 1)- Qn~^ can be analytically continued onto the lower half-plane by using the 
relation for hypergeometric functions (15.3.7) in |ll|, which transforms into the following 
relationship between q^ and q4 

q ( n +) (r; C) = C +1 qlT\r; C) + 2«p(r; ()p„(r; C) (A15) 
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Then the formulae for the analytical continuations of gi^i 2 an d gn^n 2 are 



TT ft — 1 \ f) ni 

9^l(r; 7) = 9i{;l(r; 7) - - \^-J ^ p(r; C)ftn(r; CK 2 (t; C), (A16) 
7) = <?ii(r; 7) - ^ (^) ^p(g A*(t; C)Pn 2 (r; 0- (A17) 



2. Two-dimensional operators 



Note first that the completeness relations for eigenfunctions of the one-dimensional op 



erators 



()£ + 2kt + //C£ and i) v + 2kt + /iCr? may take the form (see, e 



• g., Q; 



± 



7 



1 



(A18) 



The two-dimensional operator 



is also treated analytically 



t)z + l)v + 2kt o + AtC(£ + V) 
within the context of the basis {0„, m (£, ^)}^° m=0 (iToTK In particular, in view of (IA180 . 
the inverse of the infinite matrix 



[h 5 ®I TI + l£<g)h v + 2kt l^ ®I V + fiC (Q € ® I„ + I ? <g> Q,,) 

can be represented in the form of a convolution integral 

00 

G (±) (t ; £) = ±^ / drg^\r; 7) ® g" (±) (r - r ; 7) 



(A19) 



(A20) 



where ro = - to- Finally we noted that the matrix representation (1A20I) of the two- 

dimensional Green's function operator is symmetric in £ and 77 unlike the corresponding 
formula obtained in |2J. 
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Figure 1: The paths of integration on the physical energy sheet. The gray line is the branch 
cut along the positive real axis. The poles of the integrand in (|22H which occur at Ei = — 2 ^ 2 , 
t = 1,2, ... for t < are shown as solid circles. 
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Figure 3: The path of integration C3. The solid line is the part of C3 which remains on the physical 
sheet. The part of C3 which moves onto the unphysical sheet is depicted by the dashed line. 



18 



Table I: The nu 
element 


merical results obtained for integrals along the contours C\ and C3 and for the 

of the matrix product \)& + \ 
0,0 - 




contour C\ 


contour C 3 


£\ = Wo + xi, Si = iy + xi, 
£ 3 = k 2 - iy - \ (xi + x 2 ) 


£1 = ^+iyi, £ 2 = ^ +^2, 
£ 3 = ¥-§(yi+y 2 ) 


y = 50 , = 
y = 100, 4 1} = 
y = 500, 4 1} = 


0.99996555 - i 2.7568672 x 10~ 6 
0.99996555 + ^4.8873452 x lO" 6 
0.99996539 + ^4.9623285 x 10~ 6 

yo = 100 


vf ] = 1.0000132 - % 1.0613277 x 10~ 4 



Wl = -0.50009660 - i 0.75786560 x 10~ 4 
w 2 = -0.50011010 - i 0.77244545 x 10~ 4 
1.00013285 -i 1.3211843 x 10~ 4 



w 3 

1 = 0.79996314 x 10 



Wl = 1.49998186 - i 0.36607003 x 10~ 3 
w 2 = 1.49996754 - i 0.36818990 x 10~ 3 
w 3 = 0.99998324 - i 0.42641600 x 10~ 3 



-9 



i 0.62352239 x 10' 



(+) 



0.99999770 - i 2.3034284 x 10~ 5 



0,0 
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